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We consider U(N) and SU(N) gauge theory on the sphere. We express the 
problem in terms of a matrix element of N free fermions on a circle. This allows 
us to find an alternative way to show Witten's result that the partition function 
is a sum over classical saddle points. We then show how the phase transition of 
Douglas and Kazakov occurs from this point of view. By generalizing the work 
of Douglas and Kazakov, we find other "stringy" solutions for the U(N) case in 
the large N limit. Each solution is described by a net U(l) charge. We derive a 
relation for the maximum charge for a given area and we also describe the critical 
behavior for these new solutions. Finally, we describe solutions for lattice SU(N) 
which are in a sense dual to the continuum U (N) solutions. 
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1. Introduction 



In their classic work, Gross and Witten[l], and independently Wadia[2], showed 
that lattice QCD in two dimensions contains a third order phase transition in 
the large N limit. Basically, what happens is that the partition function for the 
plaquettes reduces to a product of partition functions for the individual plaquettes. 
The theory on each one is equivalent to a d = unitary matrix model, with the 
potential given by 

V = -^ti{U + U^). (1.1) 

This theory is equivalent to zero-dimensional field theory of N eigenvalues with 
values that lie on a circle and a potential 

V = -(Wa 2 ) + e -*) - log [sin ' . (1.2) 

i i<j 

The theory becomes critical by tuning the lattice size such that the highest fermion 
of the Dyson gas lies at the top of the potential. For fixed g, as the lattice size is de- 
creased, the potential becomes very deep and the fermions move further away from 
the critical configuration. Thus in the limit of vanishingly small a, the fermions 
should only see a quadratic potential. 

Very recently there has been a revived interest in large N QCD 2 . This program 
was started by Gross [3], who was looking for a string formulation of the problem. 
Further work in this direction was done by one of the authors [4] and finally, the 
complete theory was worked out by Gross and Taylor, who showed that QCD 2 
is essentially a string theory for a target space with arbitrary genus. The basic 
insight in handling this problem is to equate U (N) or SU (N) representations with 
sums of maps of world-sheets into the target space. It turns out that the number of 
boxes for a Young tableau of a representation is equal to the number of coverings 
of the surface. For surfaces with genus g > 1, organizing the partition function 
by the number of covers also nicely organizes the sum into powers of 1/N. Hence 
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it is very easy to compute the leading order contribution in the large N limit 
for these surfaces. For the torus, there can be leading order contributions from 
any number of coverings, but it is easy to sum their contribution. However, in 
the case of the sphere finding the leading order contribution in 1/N is not so easy, 
namely because all SU (N) or U (N) representations contribute to the leading order 
behavior, and computing these contributions involves computing 1/N corrections 
to the dimensions of these representations. 

Recently, Douglas and Kazakov (DK) discovered a clever way to solve the 
problem of the sphere by treating the rows of the young tableau and the number of 
boxes in each row as continuous variables[7]. Doing this they were able to compute 
the leading behavior which led them to a surprising result: the theory contains 
a third order phase transition, similar to the case of the lattice version of Gross, 
Witten and Wadia, but for the continuum limit of the theory. 

However, the U (N) case has an interesting feature — the sum over representa- 
tions is not asymptotic. In order to correct for this, one must sum over an infinite 
number of charge sectors. While this will eventually lead to overcounting for finite 
N, the answer will at least be asymptotic in the large N limit. These different 
charge sectors can be thought of as being extra solutions to the large N equations 
of motion. These different sectors are conjugate to the U(l) instantons. Because 
of this, one should be able to find the such solutions using the DK analysis. In this 
paper we do precisely this by generalizing an ansatz described by DK. 

There has also been some recent work by Witten in QCD 2 , although from a 
different perspective [8]. He has shown that the partition function on any Riemann 
surface is given by a sum over the saddle points. An interesting question is how 
the DK phase transition appears from this point of view. 

In section two we review recent work on 2d QCD and construct a free fermion 
picture for this theory. In section three we present an alternative derivation of 
Witten's result that the partition function is a sum over saddle points. We then 
show how the continuum phase transition occurs in this dual picture of QCD 2 . In 
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section 4 we consider the new solutions for U(N) QCD 2 and generalize the analysis 
of DK for these solutions. We consider the two cases where the area is near its 
critical value, and when the area is very large. In section 5 we compare these 
solutions to corresponding solutions for lattice SU(N) QCD 2 . In section 6 we 
present our conclusions. We include an appendix with some useful equations for 
elliptic integrals. 



Let us first review how QCD 2 on a cylinder is the same as a theory of free 
fermions by showing that it can be reduced to a one-dimensional unitary matrix 
model[9-ll]. In the gauge Aq = 0, the Hamiltonian is given as 



o o 

with the overdot denoting a time derivative. The Aq equation of motion is now 
the constraint 



2. Review of QCD 2 



L 



L 




(2.1) 



D l F lQ = d 1 A 1 +ig\A 1 ,A 1 ] = U. 



(2.2) 



Define a new variable V(x), 



V(x) 



L 



(2.3) 



where 




p e igf a dxA K 



(2.4) 



Then (2.1) can be written as 



djV(x) = 0, 



(2.5) 



so V(x) is a constant. Thus V(0) = V(L), which implies that 



[W,A 1 (0)] = 0, 



(2.6) 



where W = W and we have used the periodicity of A\ in x. 
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From the definitions (2.3) and (2.4), we find the relation 

L L 

W = ig J dxW^A x {x)W^ = ig J dxV(x), (2.7) 



and therefore using (2.5) and (2.6), we derive 

W = igLWAi(0) = igLAi{<d)W. (2.8) 

(2.8) then implies that 

[W, W] = 0. (2.9) 
Because V(x) = V(0), A\(x) satisfies 

= WgA^wS. (2.10) 

Thus, using this relation along with (2.8), we can rewrite the Hamiltonian in (2.1) 

as 

H = ^(W^W) 2 . (2.11) 

If the gauge group is U(N), with the U(l) coupling given by g/N, then (2.11) is 
the Hamiltonian for the one-dimensional unitary matrix model. The constraint in 

(2.9) reduces the space of states to singlets[15]. Hence, the problem is reducible to 
the eigenvalues of W. 

Upon quantization, this problem is equivalent to a system of iV nonrelativistic 
fermions living on a circle, with the Hamiltonian given by 

0<9 l <2n. (2.12) 

' i=l * 

The fermionization is due to the Jacobian of the change of variables from W to its 
eigenvalues, introducing the Vandermonde-type determinant in the wavefunction 



■5 



of the states, which in the unitary matrix case reads 



sm 



Qi — 6j 



(2.13) 



i<j 



2 



J|. e t(AT-l)0i ' 



where A(Aj) = rii<j(^ — -\?) * s the standard Vandermonde determinant. Notice 
that each factor in (2.13) is antiperiodic on the circle. Thus, if N is even the 
fermions have antiperiodic boundary conditions. Likewise, if N is odd they have 
periodic boundary conditions. This can be understood in terms of transporting 
a fermion once around the circle, passing by iV — 1 other fermions along the way 
and therefore picking up iV — 1 minus signs. Hence, in either case, the ground 
state is built by filling all states with wave numbers between —N/2 + 1/2 and 
N/2 — 1/2, inclusive. Subtracting off the ground state energy, one easily sees that 
this spectrum reproduces that found for the different representations of U(N). 

If the gauge group is SU(N), because A\ is now traceless W will also obey the 
condition det W = 1. Therefore the center of mass coordinate for the fermions is 
absent and we must mod it out of the theory. This means that we need to identify 
states in which all fermions have their momentum shifted by the same amount. 
Moreover, we must subtract the energy of the center of mass from the energy of 
each state in the theory. 

Now consider the partition function for the sphere. In terms of the fermions, 
we want a matrix element that corresponds to the sphere topology. The obvious 
thing to do is to map the end points of a cylinder to two points and therefore 
W = 1 there. In terms of the fermions, this corresponds to computing the inner 
product of N fermions whose position is at the point x = on the circle at time 
t — with the fermions at the same point at time t = T. This inner product Z is 
given by 



Z= (xi-t = T\ yi -t = 0), 



(2.14) 



in the limit that all Xj 



and yi — > 0. 
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To calculate Z, one can insert a complete set of momentum states for each 
particle. The wave function at t — and t = T must be antisymmetric under the 
exchange xi — > Xj or m — > yj. Therefore, the matrix element in (2.14) has the 
factor 

det|e ipi ^|det|e _ipi »|. (2.15) 
As Xi — > and yi — > 0, the factor in (2.15) approaches 

Y[(xi - Xj)(yi - yj)(pi - pj) 2 , 
i<j 

The matrix element in (2.14) is therefore given by 

Z = C 'EII (pi ~ p ^ Xi ~ Xj){yi -yj) e M-b 2 LTY,ph, (2-17) 

Pi i>j i 

where C is an unimportant constant. Not surprisingly, this term approaches zero 
in this limit. But in order to find the sphere contribution, one should notice that 
the fermion wavefunction at the end points are more singular than a 5-function, 
namely 

il>(xi) = A( Xi )S(W - 1) = -r^S(xi) (2.18) 

where a factor of 1/A 2 (xi) was produced by the change of variables from W to 
Xi. Therefore, it is necessary to divide the expression in (2.17) by these extra 
Vandermonde determinants, that is, 

Y[(xi-Xj)(yi-yj), 

leaving a finite expression. Hence the sphere partition function is 

Sphere = C ^ \[{ Pi - Pj f exp(-^ 2 LT (2.19) 

Pi i>j i 

We can compare this to the sphere partition function of Migdal and Rusakov[12,13],| 



(2.16) 
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which is given by 



Zmr = ^2{d R ) 2 e^{-9 2 AC 2R ) 



(2.20) 



R 



where the sum is over all representations of U(N) or SU (N) and cIr is the dimension 
of the representation and C 2 r is the quadratic casimir. The correspondence of the 
fermion states with the U(N) represenatations is as follows [10]: If we describe a 
representation by a Young tableau, then the number of boxes in row i, ni is the 
momentum shift of the fermion with the i th highest momentum above its ground 
state value. In terms of boxes, the casimir is given by 



which one can easily checked is reproduced by the fermions after subtracting off 
the ground state energy. The dimension of the representation is given by 



The first product in the second line of (2.22) is a representation independent term 
and is thus an unimportant constant. The second term is just pj — pi for the 
fermions. The total momentum is the U(l) charge for a representation. Hence we 
find full agreement with the result of Migdal and Rusakov. 




(2.21) 




(2.22) 



i>j i>j 
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3. Classical Solutions 



After a cursory inspection of (2.19) it would appear that Sphere is simply 
the partition function for a d = matrix model in a quadratic potential. Unlike 
the lattice case, the potential never turns over, so one might not expect a phase 
transition. As was shown by Douglas and Kazakov, this is not correct. The point 
is that the variables pi that appear in (2.19) are discrete, hence the density of 
eigenvalues will be bounded. When this bound is reached, a phase transition occurs. 
Thus, in the strong coupling phase we simply have condensation of the fermions 
in their momentum lattice, which gives a very simple physical understanding of 
the phase transition mechanism. The critical value of the area is reached when 
the density of pi, as given by the Wigner semicircle law which is valid in the 
continuoum, reaches somewhere the lattice bound, namely one, thus reproducing 
the DK result. 

There should be a classical field configuration, termed master field, which dom- 
inates the path-integral in each phase. In fact, as was shown by Witten[8] using 
the localization theorem of Duistermaat and Heckman (DH), the full QCD 2 path 
integral can be written as a suitable sum over classical saddle points. In what fol- 
lows we will give a very simple demonstration of Witten's result using the fermion 
picture and identify the classical configuration which dominates, that is, the master 
field. 

The key observation is that the exact propagator of a free particle is propor- 
tional to the exponential of the action corresponding to the classical (straight) path 
connecting the initial and final points. Since QCD 2 on the torus is equivalent to N 
free fermions, its partition function will also be given by an appropriate classical 
path of the free fermions. The things to be taken into account, however, are 

i) The particles live on a circle; therefore there are several possible classical 
paths for each of them, differing by their winding around the circle with fixed 
initial and final positions. 
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ii) The particles are fermions; thus one should also consider paths where the 
final positions of the particles have been permuted, weighted by a fermionic factor 
(— ) c , where C is the number of times the paths of the particles cross. 

The total partition function will then be the (weighted) sum of the actions of all 
these classical configurations. Since to each path corresponds a (diagonal) matrix 
W(t), and to that (up to gauge transformations) a classical field configuration 
satisfying the field equations of motion, this is the sum over saddle points of the 
action of Witten. 

For the sphere the same picture holds, with the difference that all paths start 
and end at the point x — 0, and that each path is further weighted by an extra 
factor, due to the division by the Vandermonde determinants as explained in the 
previous section. This is, again, the sum over saddle points of Witten, and the 
extra weighting factors are the determinants which appear in the DH theorem. 
These paths are characterized by their winding numbers {rii} (up to permutation) 
and thus this is a sum of the form 



where w{rii) are the (as yet undetermined) weighting factors. 

The easiest way to obtain the full expression in (3.1) is to Poisson resum (2.19). 
Using the formula 




(3.1) 




(3.2) 



n 



n 



where f(x) is any function and / is its Fourier transform, we obtain 




(3.3) 



where 
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To find the Fourier transform appearing in (3.4), we first note that 



Fi = J Hdpie-^i^Aip^eM-i^YlPi) = CA^) exp(-^- J>f). 

i i i 

(3.5) 

To prove this, notice that 

Fi = A(-d pi ) J Hdpie-^^eM-h^Pi) = exp(-^- J>?). 

i i i 

(3.6) 

polynomial of degree N(N— 1)/2; moreover it is completely antisymmet- 
ric in X{. Therefore, up to a normalization, it is the Vandermonde. The constant 
C in (3.5) can be found explicitly, it is however irrelevant for this discussion since 
it will amount to an overall coefficient in the final result. Using the convolution 
property of the Fourier transform of a product, in combination with (3.5), we find 

F 2 (xi) = (FKSFiXxi) 

= cjU dyM X -^) A(^p) exp ( - ^ $> 4 + Vl f + (xi - Vi f 

= c ex p(^4r x ^ I II d y* - 4) ex p ( - Yl y?) 

y i J i i<j y i 

(3.7) 

where and yij = yi — yj. Substituting (3.7) in (3.3) we recover 

expression (3.1), with the weights w(rii) given by 



C fX\d yi JJ(yJ _ 4) exp ( _ _J_ £ y ^ . (3.8) 



The above determines the expression of the sphere partition function in terms of 
classical saddle point configurations. 

In the large N limit one particular classical configuration in (3.1) should dom- 
inate. To determine this, let us ignore for the moment the fact that the rii are 
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discrete and replace the sum in (3.1) with an integral. Then we can bring the 
expression for Z sp h e re hito the form 

Sphere = C I J[ ds^A^) A(ij) exp ( - \N + t?)) , (3.9) 

where we rescaled g 2 to g 2 /N, in order to have a nontrivial large iV limit, and 
changed variables to s = (y + 2ixn) / ^/2g 2 LT, t — (y — 2nn) / \/2g 2 LT. We see that 
in (3.9) all explicit dependence on the coupling constant and area have disappeared. 
The area enters this picture indirectly, through the discreteness of n{. (In fact, if 
it were not for this discreteness the above integral would vanish.) Due to this, the 
plane (s, t) is discrete in the s — t direction and consists of parallel diagonal lines 
with distance D = 2n/ v /g 2 LT. If this spacing is such that the distribution of s 
and t contributing to the above integral in the large N limit is entirely within the 
lines n — 1 and n — — 1, then only the sector n = contributes. This signals a 
phase transition. Unlike the quantum case, there is no exclusion principle for the 
windings rij and thus no maximum density to be saturated. The phase transition 
in the classical expansion is more like a Bose condensation to the ground state 
n = 0. 

To estimate the critical area, substitute the Vandermondes in (3.9) with their 
absolute value, which corresponds to putting the product of differences appearing 
in (3.8) into absolute values. This has no effect for the rii = term, while it 
overestimates the contribution of the other sectors. (3.9) then becomes the product 
of two independent integrals involving a gaussian factor and one power of the 
Vandermonde. The large iV saddle point of these integrals is found in a standard 
way and the distribution of S{ and t{ is, in fact, again a Wigner semicircle with 
radius R = \[2. This on the (s, t) plane creates a "square" distribution with side 
2R in each direction. When this square lies entirely within the lines n = ±1, 
the sectors n ^ do not contribute and we are in the Boson condensate phase. 
For an estimate of the critical area, assume that the transition occurs when the 
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corners of the square start touching the lines n — ±1, that is, when the half- 
diagonal becomes D. Putting 2R\^2/2 = D we obtain g 2 A CT a = it 2 . Somewhat 
surprisingly, our estimation has given the precise result of Douglas and Kazakov. 
It is not entirely clear from this argument why the transition should occur precisely 
at this point. 

In any event, we see that the classical configuration in the weak-coupling (small 
area) phase is rii = 0. This corresponds to the master field A = 0, up to gauge 
transformations. As a check, notice that in this phase the discreteness of pi in the 
quantum expression is irrelevant (no saturation) and thus we can substitute the 
sums in (2.19) by integrals. This, upon resummation, becomes the ni = term 
in (3.1), as can be seen by the expression for w(0) obtained from (3.8). Only in 
the strong coupling phase will we have a nontrivial master field. This is implicitly 
determined by (3.1) and (3.8), although its explicit expression is not known. 



4. U(l) Sectors 

There is one difficulty with the string picture as it now stands for the U(N) 
gauge group, namely, the sum is not asymptotic with the exact answer in the limit 
1/N — > 0. This is because there exist states with finite energy, but which will never 
appear in the perturbative sum in (2.20). These are the states that correspond to 
iV fermions with momenta shifted by a constant finite amount. Such states have 
finite energy in the large N limit, but do not show up in a perturbative sum over 
surfaces since the corresponding Young tableau has at least N boxes. Moreover, 
these states are local minima, in the sense that in order to find a state with lower 
energy, it is necessary to shift the total momentum by a large amount. 

We can rectify this problem by including other sectors in the sum. For finite 
N, this will eventually lead to overcounting, but the answer will be asymptotic. 
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To this end, let us define Z, 




(4.1) 

where we have basically added m boxes to each row. After some manipulation, Z m 
reduces to 

Z m =J2 4e-^ 1+2 ™/^ n «e-w^e^, (4.2) 

Reps 

Hence we see that Z m has the same form as Zq, except that the area term that 
appears in front of n has been shifted and there is an extra factor of m 2 in the 
energy. But the n term is the same, meaning that the Gross- Taylor rules are 
exactly the same as in the Zq case, except that the area that will be used for the 
Nambu-Goto term has been renormalized, with a different shift for the chiral and 
antichiral sectors. The sum 

Z = Y J Z m (4.3) 

m 

is now an asymptotic sum for the QCD 2 string. 

We can interpret the sum in (4.3) as a sum over different classical string so- 
lutions, since around each sector m there is a sum over 1/N, the string coupling. 
These different sectors are basically the conjugates of the U(l) instanton sectors. 
We can Poisson resum (4.3), giving the expression 

Z = , /— d 2 R e ( 2 + n > n -e 2N n R e -iw^e ^a. (4.4) 

» 3 m Reps 

The m 2 term now has a factor of 1/g 2 in front of it, which is the contribution one 
would expect from U(l) instantons. Furthermore, the area that now appears in 
the Nambu-Goto action is now complex*, and there is also an additional n 2 term. 

* This might be useful for formulating QCD 2 as a topological field theory. We thank S. 
Cordes and C. Vafa for informing us of their work in this area. 
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This last term has been discussed before in the SU (N) case and was attributed to 
contributions from tubes and handles on the world-sheet. 



Since there are other string solutions corresponding to different values of m, 
we might expect to be able to generalize the work of Douglas and Kazakov to find 
the leading order contributions around these solutions. To this end, consider the 
solutions to the equations of motion for the path integral in (2.19), 



where we have absorbed the QCD coupling into the area and rescaled it by a 
factor of N. In the large N limit we can define the new variables x = i/N, and 
h(x) = Pi/N, giving the new equation, 



where u(X) is the density of eigenvalues, u(X) = dx/dX. But unlike the eigenvalues 
of a matrix model in a quadratic potential, the pi are discrete, therefore the density 
u(X) is bounded, satisfying u(X) < 1. 

DK treated this problem by dividing up the possible real values of A into three 
regions, where u(X) = 0, < u(X) < 1 and u(X) = 1. They choose the ansatz that 
the first region occurs for |A| > a, the second for b < \X\ < a, and the third for 
|A| < b. If we define u(X) to be the density of eigenvalues minus the contribution 
from the region where u(X) = 1, then (4.6) can be rewritten as 



The problem has been reduced to a two cut eigenvalue problem. To solve this, 




(4.5) 




(4.6) 




(4.7) 
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define the function 



f(h) = J dX^L. (4.8) 

f(h) then must satisfy 



where the contour surrounds the cuts from the square roots, but not the cut from 
the log nor the singularity at h. Pulling the contour back, (4.9) leads to the 
equation 

b 

A . h-b rz — — — — — r . i 



f(h) = h\ + log ^ - ^j{a 2 -h 2 ){b 2 -h 2 ) I d 



(h-s)^(a 2 -s 2 )(b 2 -s 2 )' 

(4.10) 

Expanding about large h, one can then find relations between a, b and A. 
Before doing this, let us go back and reexamine the DK ansatz. They of course 
make the reasonable assumption that the solution to the equations of motion is 
symmetric about h — 0. This is quite sensible for the ground state, but should 
not be true for solutions that correspond to sectors with nonzero values of m. 
Therefore, let us relax their ansatz somewhat and assume that < it (A) < 1 for 
the regions b < A < a and d < A < c, where a, b, c and d are to be determined. 
Then we can then proceed as before, reaching the equation 
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2 s h-c J (h - s)y/{a - s){b - s)(s - c){s - d) 
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From (4.8) and (A. 6) one finds that the density of eigenvalues is given by 



= V( a -A)(A- 6 )(A-e)(A-j) , ds 



7T J (A - s)y/(a - s)(b - s)(s - c)(s - d) 

c 

2 y/(a-X)(X-b)(X-c)(X-d) 



7r(A-c)(A-d) y/(a-c)(b-d) 

x ((c - <0n( *~°*~* q) + (A - c)K(g)) , 

(4.12) 

where K and II are the complete elliptic integrals of the first and third kind. 

Let us now expand the righthand side (4.11) in powers oil/h and compare it 
to the expansion in (4.8). Matching the terms of order h gives the equation 

b 

4 - / ds^ 1 = 0. (4.13) 

2 J y/(a-s){b-s){s-c)(s-d) 

c 

Using equation (A.l) in the appendix gives 

A = 4 -^>, (4.14) 

P 

where 

' (a — d)(b — c) 



p = y/(a-c)(b-d), 



(a — c)(b — d) 



Note that for the case c = —6, a = —d, the modulus q that appears in (4.14) is 
different than the corresponding equation in [7]. 

Matching the hP terms in the expansion of (4.8) and (4.11) leads to the equation 



b 

i f , a + b + c + d f , 

^ / ds — = = I ds 



a/(<2 — s)(b — s)(s — c)(s — d) J ^J(a — s){b — s) (s — c)(s — d) 

(4.15) 
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Using (A.l) and (A. 2), we can rewrite this equation as 



g^n(o,,) = fl + * + C -' g(,), (4.16) 
P 2p 



where a = y I n [7], the /i° equation is trivially satisfied due to symmetry. 
Matching the /i -1 term in the two equations gives 



/ 



dsu(s) + b — c = — J ds 



a/(o — s)(b — s)(s — c)(s — d) 

c 
b 

if, s(a + b + c + d) 



+ 5 / ds 



c 

b 



+ \j d s 



a/ (a — s)(b — s)(s — c)(s — d) 

a 2 + b 2 + c 2 + d 2 - 2ab - 2ac - 2ad - 2bc - 2bd - 2cd 



yj \a — s)(b — s) (s — c) (s — d) 

(4.17) 

The integral on the lefthand side of (4.17) is the total number of eigenvalues minus 
those in the region between c and b, divided by N. But the number between c and 
b is just N(b — c), thus the lefthand side is unity. Using (A.l), (A. 2) and (A. 3), 
and then using the equation in (4.16), we can reduce (4.17) to 

(a — b — c + d) 2 , . . . , 

1 = ^ L K (q) + pE(q), (4.18) 

4p 

where E(q) is the complete elliptic integral of the second kind. 
The next term in the expansion in (4.8) is 



h 2 J dsu{s)s. 



But this is just the sum of the momenta for this particular solution coming from the 
regions where the density is less than one, divided by iV 2 . This is its contribution to 
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the U(l) charge divided by N 2 . Let us call this rescaled charge Q. Using equations 
(A.1)-(A.4), (4.16) and (4.18), the h~ 2 expansion leads to the equation 

Q + Ub 2 -c 2 )=Q = ^ + ^XYZ, (4.19) 



where 



W = a + b + c + d X = a-b-c + d 

Y = a + b — c — d Z = a — b + c — d. 



(4.20) 



(b 2 — c 2 )/2 is the contribution to the rescaled charge from the eigenvalues that sit 
between c and b, hence the lefthand side of (4.19) is the total rescaled charge. The 
charge for the lowest energy state in sector m is mN , hence Q = m/N. 

Finally, we can find the specific heat, which is the next term in the expansion of 
f(h). Using (A.1)-(A.5), (4.16) and (4.18), we find, after a fair amount of algebraic 
manipulation 



)s 2 



F'(A,Q) - 1/24 = J dsu(s) 

-{3W' Z + X 1 + Y' 1 + Z l ) 



48 



+ ^M(3py2 x 2 + q WX YZ + Y 2 Z 2 + 2X 2 Y 2 + 2X 2 Z 2 + X 4 ) 



Yip 



(QWXYZ + Y Z + X l Y l + X l Z l ) 



+ ^-(3W 2 + X 2 + Y 2 + Z 2 ). 
48 

(4.21) 

Given the area A and the charge Q, in principle, one should be able to solve for 
a, b, c and d using the four equations (4.14), (4.16), (4.18) and (4.19). However, 
we should not expect solutions for all possible values of A and Q. For instance, 
from DK, we know that there are no solutions for Q = and A < ir 2 . This of 
course is the weak coupling regime. But for nonzero values of Q, the minimum 
value of A should be higher. In fact, in the infinite area limit, the upper value of Q 
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is bounded. Naively, this should happen for the sector when the fermion with the 
smallest or largest momentum is zero. We will see that the naive answer is correct. 

To further analyze the problem, we will consider two regions. The first corre- 
sponds to the area near its critical value for small values of Q. The second region 
occurs for very large values of the area. Let us first consider what happens at the 
critical point. In this case b — c, thus a = q = and U = K = E = it/2. From 
equation (4.16) we find that a = —d at the critical point. To determine b, let us 
go back to equation (4.12). Since b = c, we find that u(X) is given by 

«(A) = J°~l J—^ [(A - b) + (b + a)} 
V a + A Va 2 - b 2 

[ a 2 - A 2 
~ V a 2 - b 2 ■ 

If b 7^ 0, then for A 2 < b 2 , u(X) > 1. But this violates the ansatz that the density 
is less than or equal to unity. Therefore, we must choose 6 = 0. Using (4.18) and 
(4.19) we then see that a = | and that the charge for this solution is zero. 

Now consider small, but nonzero values for Q. In this case, b — c must be 
nonzero. To this end, let 

e = b - c, 5 = b + c (4.23) 

Using (4.18) and the asymptotic expansions in the appendix, we find the leading 
correction to a from its critical value is 

Aa = -^(e 2 + 25 2 ). (4.24) 

Since this correction is of order e 2 and not e, it will not contribute to a + d in leading 
order. This leading order correction can be found from (4.16) and the asymptotic 
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(4.22) 



expansions in the appendix. A little algebra shows that 



a + d= ^e 2 S. (4.25) 
We can now use (4.19), (4.23) and (4.25) to find Q. Let us rewrite (4.19) as 

Q = (a + d) (^ + ^-(a-d+b-c)(a-d-b + c)j 

(IK \ (426) 

+ (b + c) I - - —{a - d + b- c)(a- d-b + c) J . 
V4 16p J 

The factor multiplying a + d is to leading order 1/2. The term multiplying b + c 
is actually much smaller. In fact, this term is third order in e and is given by 
— 7r 3 e 3 /1024. Hence, the leading contribution to the charge actually comes from 
the a + d term and is 

Q = ^e 2 5. (4.27) 

The charge that appears in (4.27) is limited by the maximum value of 5 given 
e. This bound is determined by enforcing the ansatz that u(X) < 1. The region 
where this ansatz might be violated is where A ~ b or A ~ c. Let us consider the 
case where A is near b. Examining equation (4.12), we see that the first modulus 
in the elliptic integral of the third kind approaches unity as A approaches b from 
above. Therefore, if we substitute the asymptotic expansion for II in (A. 11) in 
(4.12), we find that the density is given by 



, u 2 1 (a-b)(X-b) \ n , . _ Aa-c){b-d) 

u(X) = -^J l (a _ c ^_^ (6 " c)K {q) + (c - d)K - 



+ (c - d)\ " jyff " - C j |+Q(A- 6) 



a-6)(c-d) V (A-6)(c-d) 



(4.28) 
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Hence, in order to ensure that the ansatz is satisfied, it is necessary that the relation 



°—^E{q) - K(q) > 0, (4.29) 



be upheld. Using the asymptotic expansions and the values for b, c, a and d given 
by (4.23) and (4.24), we find that 

^ b E( q )-K(q)^^(25e + e 2 ). (4.30) 

Therefore, in order for (4.29) to be satisfied, we must have 5 > — e/2. We can also 
derive the constraint that 5 < e/2 by examining u(X) as A — > c. Therefore, we 
must satisfy 

\5\ < e/2. (4.31) 

We can rewrite this constraint in terms of the charge and the area. From (4.14) 
and the asymptotic expansion, we have the relation 

3tt 4 

A-A c = — (e 2 + 25 2 ), (4.32) 

where A c is the critical value for the area. Hence, using (4.27), (4.31) and (4.32), 
we have that the maximum allowed charge sector for a given area near its critical 
value is 

Qma X = ^-(A-A c f 2 . (4.33) 

This relation is a little reminiscent of the maximum charge of a Reissner-Nordstrom 
black hole. 

Now we wish to examine the behavior deep in the strong coupling regime, which 
corresponds to large values of the area. In this case, we expect b to approach a and 
c to approach d. The values of a and d are determined by the charge of the sector 
that we are considering. For this behavior q — > 1, and thus q' — > 0. At this point it 
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is convenient to rewrite U(a, q) in terms of elliptic integrals of the first and second 
kind. Using the relation in the appendix, equation (4.16) can be written as 

-E(q)F(6, q) + K{q)E{6, q) = - + & ~ - + - K(q), (4.34) 

2p 

where sin^ = Using (4.18) and the asymptotic expansions in (A. 12) and 
(A. 13), we find that in this region 

a - d « 1 (4.35) 

which is expected since the the integral of -u(A) should be unity and for almost all 
values of A between a and d, u(\) = 1. Plugging in leading order asymptotic expan- 
sions in (A.12)-(A.15) into (4.34) and invoking (4.35), we then find the following 
approximate equation 



. l + sin# (a-c)(b-d) 

-log — + logW 7 7~w -7r->logW -walog.. . 

cos V{a — b){c — d) \a — b y(a — b)(c — d) 

(4.36) 

If we let a — b = e and c — d = e M , then (4.36) gives 



1 



1 + n 



(4.37) 



Since < /x < 00, we see that the possible values of a range from to 1. Of course 
what this means is that no local solutions exist if the fermions are shifted such 
that all of the momenta are greater than or all are less than 0. This should not 
come as a big surprise, since once these limits are reached, then there are small 
deformations of the fermion momenta which lower the energy of the state. 

The charge for these solutions is dominated by the W term in (4.19), since X 
Y and Z are all small. Clearly in the limit b — > a and c — > d, the charge approaches 

Q^ a -±± = a-\. (4.38) 
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5. Correspondence with Lattice Models 



Douglas and Kazakov have remarked that the phase transition for the contin- 
uum model is similar to the phase transition that occurs for the lattice. In both 
cases the phase transition is third order and the equations of motion for the eigen- 
values are given by a two cut model. In some sense, these two situations are dual 
to each other, with the weak coupling region of the lattice model acting like the 
strong coupling region of the continuum case. 

However, it would appear that this correspondence breaks down when we con- 
sider the solutions with nonzero values of Q. There are no corresponding solutions 
for weakly coupled U(N) on the lattice. But a little more thought shows that 
correspondence is between continuum U(N) and the lattice SU(N) and vice versa. 
For instance, the continuum SU (N) case does not have these extra solutions, since 
the U(l) charge is not a degree of freedom. In terms of the fermions, the center 
of mass coordinate is modded out. Hence shifting all the fermion momenta by 
the same amount gives back the same state. On the other hand, the SU (N) case 
has an extra term in the casimir, n 2 /N 2 , where n is the number of boxes in the 
representation. But this term does not survive the scaling limit, so we can safely 
drop it. 

For SU(N) on the lattice, the eigenvalues sit in a potential described by (1.1). 
However, unlike the U(N) case, the center of mass position for these eigenvalues 
must satisfy the constraint that 



where m is an integer. Clearly, the U(N) classical solution is the same as the 
lowest energy state for the SU(N) case, which will have m — 0. But suppose we 
consider a case where < m < N/2. Since m is an integer, we cannot smoothly 
deform this to the m = classical solution. Hence each value of m must have 
a classical minimum. We could also have values of m which are less than 0, but 




(5.1) 
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shifting the total position by 2ttN just shifts all eigenvalues around the circle, 
hence these correspond to the same solution. The maximum value of m for a given 
weak coupling is determined by the value which puts an eigenvalue at the top of 
the potential. Clearly as the coupling becomes stronger, the number of classical 
solutions will decrease. 

Let us discuss this situation in a little more detail. We can impose the con- 
straint in (5.1) by inserting a 5-function into the path integral with the form 

*(5> - 2*™) = ^-Le-sGW. (5.2) 

Hence, using (1.2) we find the equation of motion 

N(3 sm9 t - cot ^ o 6j ~-(J2 9i ~ 2?rm ) = °- ( 5 - 3 ) 

where (3 = l/(a 2 g 2 ). Under the usual rescaling, one ends up with the equation 

psin9-\J n{<f))(f> - 2nm/N) = P J ^(0) cot (5.4) 

where fj,((f>) is the density of eigenvalues. As e — > 0, the last term on the lhs of (5.4) 
is divergent unless the integral is very close to 2-kitl. To this end, let us replace 
this entire term by a constant a. Next define the function f{6), 

f{9) = f dM<P) cot 9 -^. (5.5) 

Given this definition, the function is analytic everywhere in the strip —n < re6 < ir, 
except for a cut along the real line. f{6) is clearly invariant under 6 — > 9 + 2n. 

f(9) can be solved for by analyzing its behavior as Im# — > oo. In particular, 
using (5.5) the large Im^ behavior for f(9) is 

f{9) = ~ij d<j>n(<i>) - ie ie J #e _i ^(0) + 0(e 2id ). (5.6) 
Let us assume that the eigenvalues sit on the strip A — a < \ < A + a. Then from 
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(5.4) and (5.5), /(A ± ie) satisfies 



f(\±ie) =/3sin(A) + a ^ 2ivr/i(A). (5.7) 
Choosing the periodic function 



9 + A I . 2 9 - A . 2 a 



/(0) = (3sm9 + a-2(3cos——\sm z — sin 2 -, (5.8) 

Zi V — — 

one finds that the asymptotic behavior is matched if 

/3 sin 2 ^ cos A = 1, (5.9) 



and 



/3cos 2 ^ sin A = -a. (5.10) 



From (5.7) and (5.8) one learns that the density of eigenvalues is given by 



P _ X + A /_. 2 a . 2 A - A 



/i(A) = -cos^— y S in 2 --sin 2 ^— . (5.11) 

In order that the density of states is positive, one must satisfy the inequalities 
A + a < n and A — a > —n. One can also easily show that the integral of sin0 
weighted by the density of eigenvalues satisfies 



/ 



/j,(4>) sin0 = cos 2 - sin A. (5-12) 
2 



Hence A basically measures the anisotropy of the solution about the bottom of the 
well. 

As (3 — > oo, it is clear that the allowed solutions for A lie anywhere between 
— 7r/2 and it/2. This is similar to the situation in the previous section, where any 
charge between —N/2 and N/2 is allowed as the area approaches infinity. However, 
near the critical point, only a small range for A is allowed. Clearly, from (5.9), the 
critical value of f3 is f3 c — 1. Therefore, the only allowed A at the critical point is 
A = 0. 



2G 



Moving slightly away from the critical point, one has a = n — e and a maximum 
value of A given by A = e. At the maximum, the integral of the angles weighted 
by the density of eigenvalues is then given by 



A+a 



AM / J.\ J. P f jmm <^ + A /• 2« • 20~ A 

a<f>fj,{<p)<f> = — / ct00 cos — - — ysm — — sin — - — 



A-a 



(5.13) 



P ■ a f i i , ■ i I ■ 2 a ■ 2 

sin A / a00sm0ysm — — sin — , 



where we have shifted the integration variables and have used the fact that the 
integral over the density of states is one. For small e we find that the leading order 
term in (5.13) is 

/ rf0/i(0)0 = ie 2 A log - + 0(e 2 A). (5.14) 
J 8 e 

f3 is the inverse area, hence we find that the area for the plaquette satisfies 

A c -A = i(e 2 + 2A 2 ). (5.15) 
Therefore, the maximum charge in terms of the area, behaves like 

Qmax ~ (A c - Afl 2 log -^-j. (5.16) 

Hence, unlike the continuum case, the maximum charge in terms of the area has 
a scaling violating piece. However, the integral over sin0 has a maximum that 
behaves like (A c — A) 3 / 2 , which is closer to the behavior of the preceeding section. 
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6. Discussion 



In this paper we have given an alternative derivation of Witten's proof that 
the QCD partition function is given by a sum over saddle points. We then showed 
how the sphere phase transition occurs in the dual picture of QCD string theory. 

We next demonstrated how to find new solutions for the U(N) QCD 2 sphere 
partition function corresponding to the different U(l) sectors. We have also shown 
how these solutions relate to solutions found for the lattice versions of SU(N). In 
some sense, it is surprising that the nonzero charge solutions can actually be found. 
One might have expected that since there is nothing in the machinary of section 
4 that explicitly states that the eigenvalues lie on the lattice, then nothing should 
prevent the nonzero charge solutions from sliding down to the absolute minimum. 
But actually, this information is in there, because of the log term that appears in 
(4.9). This term states that there is a region where the eigenvalues have a constant 
density, that is, they lie on the lattice. The boundaries of this region essentially 
add another degree of freedom. This then allows us to find solutions with nonzero 
charges. 

It is hoped that the results presented here might have some use in investigating 
random matrix model theory. Perhaps one can consider cases where the entries of 
the Hamiltonian are restricted to be integers. One might then find similar behavior 
to that shown here. 



In this appendix we present some useful formulae of elliptic integrals. These 
are taken from or are easily derived from formulae in [16]. The first such equations 
are 



APPENDIX 



b 




(A.l) 



c 
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b 

[ ds ; S = -((c- d)U(a, q) + dK(qj) , (A.2) 

J J(a-s)(b-s)(s-c)(s-d) ' V ' V V V ' 

c 



c 



a/(<2 — s)(6 — s)(s — c)(s — rf) 

i((-ad-6c + d(a + 6 + c + d))if(g) - p 2 E(q) + (a + 6 + c + d)(c - d)U(a, q) 
p\ 

(A.3) 

b 



I 



s 3 



c 



rfs 

y (a — s)(6 — s)(s — c)(s — d) 



i(2d 3 - i(c-d) (6- d) (a + 36 + 3c + 5d))fT(g) - ^(a + 6 + c + d)pE(q) + 
a + b + c + d)(a + b + c - 3d) 



P 
1 /3 



PV4 

+ (2a + 2b + 2c + 3d)d -ab-ac- bcj (c - d)n(a, q), 

(A.4) 

and 



ds- 



c 



a/(<2 — s)(6 — s)(s — c)(s — d) 



~p ~ ^ c ~ d ^ b ~ d "> t~ 5fl2 - 15&2 - 15c2 - 33rf2 

- 4a6 - Aac - 6ad - 146c - 246d - 24cd]) K(q) 

- — (aha 2 + 156 2 + 15c 2 + 15d 2 + 14a6 + 146c + 14ac + 146c + 146d + 14cd)p£(g) 

+ — ( 5a 3 + 56 3 + 5c 3 + 5d 3 + 3a 2 6 + 3a 2 cx + 3a6 2 + 36 2 c + 36c 2 + 3a 2 d 
8p V 

+ 3ad 2 + 36d 2 + 3cd 2 + 2a6c + 2a6d + 2acd + 26cd) (c - d)U(a, g), 

(A.5) 
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where 



-E(g) and n(o;,g) are the complete elliptic integrals of the first, second and 
third kind respectively. 

Another formula used in the text is 



1 



(A — s) y/(a — s) (b — s)(s — c) (s — c) 
2 



(A - c)(A - d)y/(a-c)(b-d) 



c-d)U( b b _ C d X x _ d c , q ) + (h-c)K(q) 



(A.7) 

We also use the relation between the complete integral of third kind and in- 
complete integrals of the first and second kind, 

q) = ^^-K(q) - E(q)F(e, q) + K(q)E(6, q), (A.8) 
P P 

where sin 2 9 = and F(6, q) and E(6, q) are the incomplete integrals of the first 
and second kind. 

The following asymptotic expansions are also used: 



E (q) = 1 ( 1 _ i 5 2 _ 1 4 _ J^6 + ) (A 1Q) 

w 2 I 4 y 64 y 256^ ' ' 



n(a, g) =- (\ + i(a 2 + g 2 /2) + ||(a 4 + a 2 g 2 /2 + 3g 4 /8) 

+ ^(a 6 + «V/2 + 3a 2 g 4 /8 + 5g 6 /16) + ...), 



(A.ll) 
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K(q) = log(4/g') - 2[log(4/g') + Wf + (A.12) 

E(q) = 1 + ^pog(W) " V2](g) 2 + (A.13) 

F(0, g) = log + \ [log - sin 9 sec 2 £](g') 2 + .... (A.14) 

cos# 4 cos# 

E(6,q) = sm6 + 2[log 1 + S1 " 1 - - sing] (q'f + .... (A.15) 

cos 6* 



where g' satisfies </ = a/1 — g 2 . 

Finally, there is a useful expansion for U(a, q) if at is close to unity. This is 

n (a , 9 ) = ^)-g + ^^'» + o( 1 - a »). (A,a) 
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Note added: As this paper was being typed we received a paper[17] which 
discusses some of the issues in section 3. 
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